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VARIABILITY IN THE THERMAL EMISSION FROM ACCRETING NEUTRON STAR TRANSIENTS 
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ABSTRACT 

The composition of the outer 100 m of a neutron star sets the heat flux that flows outwards from the 
core. For an accreting neutron star in an X-ray transient, the thermal quiescent flux depends sensitively 
on the amount of hydrogen and helium remaining on the surface after an accretion outburst and on the 
composition of the underlying ashes of previous H/*He burning. Because H/^He has a higher thermal 
conductivity, a larger mass of H/*He implies a shallower thermal gradient through the low density 
envelope and hence a higher effective temperature for a given core temperature. The mass of residual H 
and *He varies from outburst to outburst, so the thermal quiescent flux is variable even though the core 
temperature is constant for timescales < 10^ yr. Heavy elements settle from a H/^He envelope in a few 
hours; we therefore model the quiescent envelope as two distinct layers, H/*He over heavier elements, 
and treat the mass of H/'^He as a free parameter. We find that the emergent thermal quiescent flux can 
vary by a factor of 2 to 3 between different quiescent epochs. The variation is more pronounced at lower 
interior temperatures, making systems with low quiescent luminosities and frequent outbursts, such as 
SAX J1808.4— 3658, ideal candidates from which to observe this effect. Because the ashes of H/^He 
burning are heavier than ^^Fe, their thermal conductivity is greatly reduced. This increases the inferred 
crust temperature beyond previous estimates for a given effective temperature. We survey this effect 
for different ash compositions and apply our calculations to Cen X-4, Aql X-1, and SAX J1808.4— 3658. 
In the case of Aql X-1, the inferred high interior temperature suggests that neutrino cooling contributes 
to the neutron star's thermal balance. 

Subject headings: conduction — diffusion — stars: individual (Aql X-1, Cen X-4, 
SAX J1808.4-3658)— stars: neutron— X-rays: binaries 



1. INTRODUCTION 

The orbiting X-ray observatories Chandra and XMM 
have dramatically improved our understanding of soft X- 
ray transients (SXTs): binaries containing a neutron star 
or black hole primary and having well-defined accretion 
outbursts separated by long periods of quiescence. These 
objects are typically defined as having a ratio of outburst 
flux to quiescent flux > 1000. Two puzzles are pertinent 
to this work. The first is whether the thermal component 
of the neutron star 's quiescent luminosity is powered by 
accretion or by thermal emission from the cooling core. 
This has implications for the observed contrast in the qui- 
escent luminosity between transient black holes and neu- 
tron stars (Narayan et al. 1997; Menou et al. 1999; Garcia 
et al. 2001; Kong et al. 2002). 

Brown et al. (1998) showed that compression-induced 
reactions — electron captures, neutron emissions, and py- 
cnonuclear reactions (Bisnovatyi-Kogan & Chechetkin 
1979; Sato 1979; Haensel & Zdunik 1990)— in the iimer 
crust of an accreting neutron star release enough heat 



to power a cooling luminosity of order 10^^ ergs s^^ in 
quiescence. In the absence of neutrino emission from 
the core the quiescent thermal flux is proportional to the 
mean outburst flux (Brown et al. 1998; Colpi et al. 2001). 
Motivated by the match between the expected quiescent 
luminosity and that observed from neutron star SXT's 
in quiescence, Rutledge et al. (1999, 2000) fit archival 
ROSAT and ASCA observations of Aql X-1, Cen X-4, 
4U 1608-522, and 4U 2129+47 with realistic H (or ^He) 
atmosphere spectra and found that the emission could 
be explained as thermal emission from an area of radius 
w 10 km. Further Chandra observations of Cen X-4 (Rut- 
ledge et al. 2001b), Aql X-1 (Rutledge et al. 2001a), and 
KS 1731-260 (Wijnands et al. 2001; Rutledge et al. 2001c), 
as well as quiescent transient neutron star identifications 
in uj Cen (Rutledge et al. 2002b), NGC 6440 (Pooley et al. 
2002), and 47 Tuc X5 and X7 (Heinke et al. 2002) confirm 
that the quiescent spectra of neutron star transients is con- 
sistent with being thermal emission (effective tempera- 
ture fcBTeff > 100 eV) from a H photosphere, plus, in most 
cases, an additional hard power-law component. The ori- 
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gin of the power-law tail remains uncertain (for a review 
of proposed mechanisms, see Campana et al. 1998; also 
see Menou & McClintock 2001). 

The second puzzle is the source of the observed vari- 
ability in the quiescent emission on timescales > 1 d. In- 
deed, it was in part because of an apparent increase in the 
quiescent intensity of Cen X-4 between 1980 {EINSTEIN) 
and 1984 (EXOSAT) that led van Paradijs et al. (1987) to 
discoimt the possibility that the observed emission was 
intrinsic to the neutron star, i.e., not powered by accre- 
tion. ROSAT/HEl observations of Cen X-4 (Campana 
et al. 1997) revealed that the intensity decreased by a fac- 
tor of w 3 over 4 days. Similarly, there was a fractional de- 
crease of 40 % in the observed intensity between an ASCA 
observation and one 5 yr later with Chandra, although 
this could be attributed to variability in the power-law 
component only (Rutledge et al. 2001b). A compari- 
son of Aql X-1 observations taken with Chandra/ ACIS- 
S, ROSAT/PSPC, and ASCA (Rutledge et al. 2001a) in- 
dicates variability by a factor of 2 over a timescale of 
roughly 8 yr. Note that in this case there were several 
intervening outbursts between the different observations. 
In both cases, there was no short-timescale (< 10^ s) vari- 
ability detected (Rutledge et al. 2001b,a). More recently, 
Rutledge et al. (2002a) used the Chandra/ ACIS-S to take 
four "snapshots" of Aql X-1 after a recent outburst. The 
intensity was observed to decrease by a factor w 0.5 over 
3 months and then increase by a factor w 1.4 over 1 
month. In addition, short-timescale variability was found 
in the last observation. 

The standard interpretation is that the observed vari- 
ability is caused by fluctuations in the quiescent accretion 
rate (van Paradijs et al. 1987; Campana et al. 1997; Brown 
et al. 1998; Rutledge et al. 2001b,a; Menou & McClintock 
2001; Dubus et al. 2001). In this manuscript, we describe a 
previously overlooked cause of variability in the intrinsic 
quiescent thermal emission: a changing envelope^ com- 
position. Even if accretion completely halts in quiescence, 
the neutron star's envelope will have a different stratifi- 
cation following each outburst. This varying composition 
can change the quiescent flux by a factor of 2 to 3 for a 
fixed crust/ core temperature. 

Previous calculations of the thermal structure of a cool- 
ing unmagnetized neutron star considered the difference 
between a purely ^^Fe envelope (Gudmundsson et al. 
1983) and one composed of light elements (H, "^He, and 
^^C) overlying ^''Fe (Potekhin et al. 1997). Consider two 
hypothetical neutron stars each with a core /crust tem- 
perature of Tb — 10* K (typical of neutron star SXT's; see 
below); one star has a pure ^^Fe envelope and the other 
has H and ^He at densities < 10^ g cm^''. As noticed by 
Potekhin et al. (1997), the large difference in opacity be- 
tween a pure iron and a light element envelope means 
that the effective temperature, Tgff, for the light element 
envelope is a factor of 1.6 hotter: Teff(^''Fe) = 1.1 x 10^ K 
and Tgff (light element) = 1.8 x 10^ K. Gudmundsson et al. 
(1983) first noticed that the thermal stratification is sensi- 
tive to the opacity in the region where the electrons are 
semi-degenerate. Coincidentally, it is in this region that 



the accreted H and *He unstably ignite. As a result the 
thermal gradient through the envelope depends on the 
mass of H and *He remaining after the previous accretion 
outburst. Our calculation thus addresses variations in the 
quiescent flux from one quiescent epoch to the next, as the 
intervening outburst changes the composition and mass 
of the outermost layers of the neutron star. Changes in the 
quiescent intensity over a timescale of months can also 
occur from differential sedimentation of ions and resid- 
ual accretion. Neither of these scenarios can explain short 
timescale (< lO'' s) variability such as just observed from 
Aql X-1 (Rutledge et al. 2002a). 

This paper first (§ 2) describes in qualitative terms 
the overall stratification of a quiescent neutron star tran- 
sient. Section 3 contains a summary of the relevant mi- 
crophysics in the calculation: the equation of state (EOS), 
diffusive sedimentation of ions, and thermal transport. In 
§ 4 we describe how changing the composition and strati- 
fication of the envelope produces variations in the surface 
effective temperature. This calculation is then applied, in 
§ 5, to Aql X-1, Cen X-4, and SAX J1808.4-3658. Impli- 
cations and directions for future study are discussed in 
§6. 

2. THE COMPOSITION AND STRATIFICATION OF 

QUIESCENT NEUTRON STAR ENVELOPES 

In the absence of accretion, the thermal structure of the 
envelope is determined, over durations much less than 
the cooling timescale of the core (i.e., < 10* yr), by the 
flux equation 

d f T Y 3 

Here k — (k^^ + K^^y^ is the reciprocal sum of the radia- 
tive opacity Kr and the conductive opacity Kc, 

(2) 

3pK ' 

with K and p being the electron thermal conductivity and 
mass density. The spatial coordinate is just the column 
depth, y — pdr — p/g by hydrostatic balance. We 
use the Newtonian form of the thermal diffusion equa- 
tion: the thickness of the envelope is much less than the 
stellar radius, so that the gravitational redshift 1 -|- z w 
[1 — 2GM/(Rc'^)]"^/^ is nearly constant across the enve- 
lope and factors from equation (1). All quantities in this 
manuscript refer to proper quantities; in particular the ef- 
fective temperature as observed far away from the star is 

T'eff.oo — TesO- + Z)"^. 

For a fixed envelope stratification, the flux equation (1) 
guarantees a one-to-one mapping between Tj, and Tgff. 
The core temperature cannot change on timescales <ti 
10* yr, so if the envelope composition were constant, then 
the basal effective temperature and luminosity would be 
unchanging from quiescent epoch to quiescent epoch. For 
an accreting neutron star the envelope composition and 
stratification are not, however, fixed. During each out- 
burst, H and *He are deposited onto the surface of the 
neutron star. After accumulation of a critical column 



^ In this paper "envelope" means the outermost layer of the neutron star where the thermal gradient is significant. This is distinct from the photosphere, 
where the emergent continuum spectrum forms. 
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i/ign ~ 10* g cm^^, the H and *He burn to heavier elements 
("ashes"), and the process then repeats. As the outburst 
wanes, there is a last episode of unstable burning (a type I 
X-ray burst). Accretion after this last t5^e I burst deposits 
a residual H/^He layer of column j/; < i/ign onto the ashes 
of previous episodes of H/*He burning. The depth of the 
light element layer is effectively unconstrained, and as a 
result Tgff can vary even though T\, is fixed. 

The composition of the ashes depends on the nature 
of the H/^He burning (for a recent review, see Bildsten 
2000); for most accretion rates, the ^He unstably ignites 
in the presence of H. The H is then consumed by the rp- 
process, a sequence of rapid proton captures onto seeds 
provided by the "^He burning (Wallace & Woosley 1981; 
Van Wormer et al. 1994; Schatz et al. 1998). Reaction net- 
work calculations, both for single-zone calculations of un- 
stable burning (Koike et al. 1999; Schatz et al. 2001) and 
for steady burning (Schatz et al. 1999), find that all of the 
hydrogen is consumed and that the reactive flow reaches 
nuclei much heavier than ^^Fe. In a recent calculation, 
Schatz et al. (2001) determined that the rp-process ends in 
a closed SnSbTe cycle; the resultant ash composition has 
a mean nuclear charge (Z) — 37 and a mean nuclear mass 
{A) — 79. Following the burst, the proton-rich elements 
quickly /3-decay to more stable species such as ^"^Ru. 

The ratio of H to *He is not well determined. AAHien 
the temperature in the crust exceeds roughly 10** K, the 
H in the outer envelope is consumed on a timescale ~ 
6 X 10^ s (XcNo/0.02) by the hot CNO cycle, Xcno being 
the mass fraction of CNO nuclei. Hydrogen is also con- 
sumed (in a thermally stable fashion) when the accretion 
rate exceeds about 2 x 10"^° M© yr^^ During the out- 
burst decay, a significant fraction of the accumulated H 
can be converted stably to *He, even after the last type 1 
X-ray burst has occurred. We include this possibility in 
OUT calciilations (§ 4). We note, however, that the photo- 
sphere would still tend to be pure H, as T^s is too cold for 
the H to be consumed on a timescale over which differen- 
tial sedimentation removes CNO nuclei from the photo- 
sphere (-^ 10 s; Bildsten et al. 1992). As a result, spec- 
tral fitting of the quiescent flux will find that the ratio 
of radius to distance remains constant even if Tgff varies. 
Changes in the photospheric abundances, as implied by, 
e.g., spectral features, more likely indicate quiescent ac- 
cretion. 

3. EQUATION OF STATE, SEDIMENTATION, AND 
THERMAL TRANSPORT 

The microphysics of the envelope enters equation (1) 
explicitly through the opacity k and implicitly through 
P(y ~ P/S^'^)- Before discussing the thermal transport 
and its effect on the thermal structure of the envelope, we 
first review the different physical regimes of the neutron 
star's envelope, the equation of state, and the validity of 
treating the envelope as being composed of distinct lay- 
ers. Figure 1 shows different regimes of the p-T plane 
for an envelope composed of a pure H layer, of column 
depth y = 10* g cm"^^ superincumbent on a ^°*Ru layer. 
The top panel illustrates conditions in the H layer; the 
bottom panel does likewise for the ^"''Ru layer. In both 
panels, the circles denote the thermal structure foimd by 



solving equation (1) for Tgff = 2.1 x 10^ K (as inferred for 
Aql X-1; top curve) and T^n = 8.8 x 10^ K (as inferred for 
Cen X-4; bottom curve). We discuss the calculations for 
these sources in § 5. 

3.1. Equation of State 

The envelope of a neutron star is composed of elec- 
trons and ions. Electrostatic interactions between elec- 
trons are negligible throughout most of the envelope for 
the temperatures of interest (see Chabrier & Potekhin 
1998), and the electrons are an ideal degenerate Fermi 
gas for = /Xc/AibT 3> 1. Here /Xg is the electron chem- 
ical potential, not including the rest mass. For ^ 1, 
/ie « (ep — wigC^), where ep = {nilc^ + ppC^)!/^ jg the elec- 
tronic Fermi energy and — {3Tr^ne)^^^fi is the Fermi 
momentum. We write the electron density tig = Yep/rriu 
where is the mean nucleon mass and delimit on Fig- 
ure 1 where the electrons are degenerate with the con- 
dition tp = 10 {dashed line). The electrons are relativistic 
where pp/nteC « [Yep/ilO^ g cm"^)]^/^ > 1. 

Where the envelope is composed of rp-process ashes, 
the total number of species (each of density Uj) is likely 
quite large. Electrostatic correlations between ions are pa- 
rameterized by 



ajhT 

y5/3 



(3) 



YeP 



10* g cm" 



1/3 



10* K 
T 



aj — (3Zy/47rMe)^''^ being the ion sphere radius for species 
/. For r > Tm, the plasma is a solid; Fm is computed by 
equating the free energies of the liquid and solid phases 
(Fig. 1, solid line). We compute the ionic free energy for 
the liquid phase, 1 < F < Fm, from the fit of Chabrier 
& Potekhin (1998) and for the solid phase from the fit 
of Farouki & Hamaguchi (1993). In a recent calcula- 
tion, Potekhin & Chabrier (2000) determined that Fm = 
175.0 ± 0.4, with a further relative imcertainty of ~ 10 % 
arising from electron screening. The calculations in this 
paper are insensitive, fortvmately, to the precise value of 
Fm- 

The ions are classical for 



h V47rp(Z)2e2 



•p,ion 



4.0 X 10^ K 



10* g cm" 



1/2 



-3 \ 1/2 



2(Z) 



» 1, (4) 



where 9 is the temperature in units of the ion plasma tem- 
perature. The ratio of melting to plasma temperatures is 

(Z) /2(A) / 10*gcm- 3^ 



'M 



= 3.4 



ip,ion 



175 



30 



{Z)J 



(5) 

and so T > Tp ion wherever the ions are composed of high- 
Z species and are in a liquid phase. When 6 < 1 {dotted 
lines, Fig. 1), quantum lattice effects become important. 
As evident from Fig. 1, in a H layer the ions are typi- 
cally weakly coupled throughout, although one should 
be careful about quantum plasma effects (for a discussion. 
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see Potekhin et al. 1997). In contrast, the ions in a ^'^'^Ru 
layer are strongly coupled wherever the electrons are de- 
generate. Quantization of phonon modes is not impor- 
tant, however, except in the very degenerate layers where 
the thermal gradient is nearly isothermal. 

3.2. Sedimentation 

The calculation in this paper presumes that the enve- 
lope is segregated into layers. An order-of-magnitude 
calculation illustrates the timescale for the envelope to be- 
come stratified. In a frame co-moving with the mean ion 
center-of-mass (CM; here at rest), the continuity equation 
for species i is 

dn- 

-J.^V-i'DVni-niWi), (6) 

where D is the interspecies diffusion coefficient and if, 
and w, are respectively the drift velocity, relative to the 
mean ion CM, and the number density of species i. This 
equation neglects terms arising from thermal diffusion, 
which are generally small in dense ionic plasmas (Steven- 
son & Salpeter 1977; Paquette et al. 1986). 

The diffusion coefficient, D, and the drift velocity, w, 
are both local quantities, i.e., they do not depend on Vp 
or VT. For a trace component (species "2") in a back- 
ground (species "1") these quantities are effectively fixed, 
and are therefore related through Einstein's relation (see, 
e.g.. Landau & Lifshitz 1987), W2 = [©/(^bT)] F2 where F2 
is the force on the trace ion and the combination D / (k^T) 
is the mobility. This relation holds where the temperature 
scale height -[d(lnT)/dr]^^ ':$> Hp = p/pg, the pressure 
scale height. Since -Hp[d(lnr)/dr] = d(lnr)/d(lnp) < 
{dhxT /d]np)s < 1, this requirement is trivially satisfied. 
The force F2 is the sum of gravity and the electric field 
needed to keep the ions from settling relative to the elec- 
trons (see, e.g., Spitzer 1962; Hameury et al. 1983), 

F2 = AjMug - ZzeE 

{(A2 - ntng, non-degenerate 
(Az-^) mng, degenerate. 

The electric field is eE — niug [Ai / (Zi + 1)] where the 
electrons are nondegenerate and eE = mug{A\ /Zi) where 
they are degenerate and contribute most of the gas pres- 
svie. 

From continuity (eq. [6]) the timescale for the trace 
ions to move a distance s is Tdr = min(s^D"^, s||zv2||^^)/ 
where the first term is the diffusion timescale and the 
second is the drift timescale. These timescales are equal 
for s — H2 — A:bT/||F2||, which is just the scale height 
for the trace ions. For mass/charge ratios much greater 
than unity and where the electrons are degenerate, H2 ^ 
Hp-, therefore over macroscopic scales s ~ Hp, the rel- 
evant timescale is the drift timescale t^,. = s||zt;2|I"^ = 
sF2{kBT/D). We can evaluate t^t once the mobility is de- 
termined. For the situation in question (y ~ 10* g cm~^, 
r ^ 10* K), the background matter exists in a liquid 
phase (1 < Fi < Fi m) (see Figure 1). We follow the 
calculation of Bildsten & Hall (2001) and estimate D 
from the Stokes-Einstein relation between the mobility 
and the drag coefficient (for a liquid sphere), Ti/k^T — 



l/(47rfl2?7)- In this expression, ^2 is the radius of a 
charge-neutral (containing Z2 electrons) sphere arotmd 
the trace particle and 77 is the viscosity of the material. 
Fits to numerical simulations of one-component plasmas 
(OCP) in the liquid regime find that the viscosity is r/ w 
(0.1 g cm-i s-i)(;o/g cm-3)(a;p,i/s-i)(fl2/cm2)(ri/10)0-3, 
where ai is the ionic spacing of the background fluid 

1 /2 

and cjp 1 — [An p{Zie / AitUu)^] is the plasma frequency 
(Donko & Nyiri 2000; Bildsten & Hall 2001). 

Using the non-relativistic degenerate electron equation 
of state to relate p and p and evaulating F2 from equa- 
tion (7), we find from the Stokes-Einstein relation the mo- 
bility, 

-«1.2xl0^sg-^p^. (8) 

For brevity we scale the surface gravity, temperature, 
and density to g — gu^O^^ cm s~^, T — TylO^ K, and p — 
pslO^ g cm^^. The value D computed from equation (8) is 
comparable (within a factor of a few) to that calculated by 
Tanaka & Ichimaru (1987) for a strongly coupled OCP un- 
der the conditions of interest. It is somewhat larger than 
the value estimated with the formalism of Chapman & 
Cowling (1952) or Burgers (1969), both of which are valid 
for weakly coupled plasmas, F <g; 1 (also see Fontaine & 
Michaud 1979; Paquette et al. 1986). 

Upon substituting equation (8) into the expression for 
W2 and using Hp as a fiducial lengthscale, one arrives at 
the stratification timescale, 

r y3.9y0.3^13 1 
T ~ 1 nS c ^1 ^2 P5 (q\ 

^dr ~ [Al YuTr (AlZi - A,Z2)\ ■ 

For ^°^Ru in *He at p^ = I, Tj = 10, and gn — 1 (appro- 
priate for Aql X-1; § 5), equation (9) implies that the ^°*Ru 
ions settle in a time of roughly 2 hr. For ^^Fe in *He, the 
timescale is roughly 7 hr. For the less dense regions of the 
envelope, the diffusive time scale is of order seconds to 
minutes, so in the absence of any circulation the envelope 
quickly stratifies after the end of the accretion outburst. 
Our calculations throughout the remainder of the paper 
assume a fiilly stratified envelope. We note that in diffu- 
sive equilibrium, the boimdary between the layers has a 
thickness --^ H2 <C Hp, which justifies our approximating 
the interface as a planar surface. 

3.3. Thermal Transport 

For the temperatures and densities in the quiescent 
neutron star envelope, the relevant opacities are Thom- 
son scattering and free-free absorption. Because the out- 
ermost layers are composed of H and "^He, the ions are 
fully ionized throughout. At typical envelope tempera- 
tures, the dominant opacity is from free-free absorption, 

AtffOCner-^/'£ZfYygff,y (10) 

3 

where the Gaimt factor for species /, gff y, contains correc- 
tions for electronic Coulomb wavefunctions, degeneracy, 
and relativistic effects. There are no fits of gs that cover 
the entire (p, T) range appropriate for this problem. The 
relation between Tb and Tea is most sensitive to k where 
the electrons are semi-degenerate, however, so we use the 
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fit from Schatz et al. (1999) that is tuned to be accurate for 
^ < 10 and moderately strong Coulomb corrections, pa- 
rameterized by -4 < In 7^ = Z^e^nie/iltfkBT) < 2. This 
fit is reasonably accurate (fractional errors ~ 10 %) when 
compared against the calculations of Itoh et al. (1985, 
1991). We calculate the Thomson scattering opacity by 
using a fit (Buchler & Yueh 1976) that reproduces the non- 
degenerate limit (Sampson 1959) and includes corrections 
for the relativistic and degenerate electronic EOS. 

As t/j increases, the degenerate electrons become more 
efficient than photons at transporting heat. The electron 
thermal conductivity is given in the relaxation-time ap- 
proximation by the Wiedemann-Franz law. 



K- 



klT 



(11) 



3 m* 

where m* — ep/c^ is the effective electron mass, and r 
is the electron thermal distribution relaxation time. In 
this approximation, the relaxation time is the reciprocal 
sum over electron-electron and electron-ion scattering re- 
laxation times, T^^ = Tg"g + ''"e in- heavy element 
layer, the large values of F (see Fig. 1) may inhibit strat- 
ification, and therefore we must consider a multi-species 
plasma. Where 9 >l,a reasonable approximation (moti- 
vated by the additivity rule in multi-ionic EOS; Potekhin 
et al. 1999) is to sum over inverse relaxation times for 
each species, t^^^^ = L/ '^e y ' where t^j are the separate in- 
verse relaxation times for electron-electron and electron- 
ion (from species ;) scattering, respectively. We calculate 
the electron-electron scattering relaxation time from the 
formalism of Urpin & Yakovlev (1980), as fit by Potekhin 
et al. (1997). The inverse electron-ion scattering relaxation 
time is 



.-1 



47re* p 



niu 



(12) 



Here v-p — p^/m* is the electron velocity evaluated at the 
Fermi surface, and A^j is the dimensionless Coulomb log- 
arithmic term that originates in the integration of the scat- 
tering rate over electron phase space. To evaluate A^ j, we 
use the fitting formula of Potekhin et al. (1999), which is 
straightforward to implement for arbitrary (Z, A). 

Where 9 <1, phonon modes begin to "freeze out," and 
the additivity rule (eq. [12]) becomes suspect. In practice 
this is not typically a concern, because scattering from 
charge fluctuations (impurity scattering) becomes more 
important than electron-phonon scattering. If the impuri- 
ties are randomly distributed, then the "structure factor" 
in the integration of the scattering integral can be set to 
unity (Potekhin, private communication; also see Itoh & 
Kohyama 1993), and Ae,imp resembles that of the liquid 
(r <C Fm) phase with a relaxation time depending on the 
rms charge difference. 



imp 



Thus the scattering differs from that in a liquid by a factor 
1^-1 = ^ (14) 

(Z)2 ~ (Z)2 ^^^^ 

With the structure factor in eq. (8) of Potekhin et al. (1999) 
set to unity, we find that the resulting Ae^imp is comparable 
to the fit given by Itoh & Kohyama (1993). 



How large might Q be? The output from a one-zone X- 
ray burst nucleosynthesis calculation (Schatz et al. 2001) 
has Q/(Z)2 = 233/372 ^ 0.17. The computation of Ae,y as- 
sumes that the separate species are arranged in a lattice; 
it is difficult to imagine how this could come about in the 
case of an accreted neutron star crust. For the conditions 
of interest in this paper {p < 10^° g cm"^, T < 10* K, and 
an rp-process ash composition), t^^^ < Y.jT^J- Impurity 
scattering is therefore not dominant, unlike the case in the 
deep crust (Brown 2000; Gnedin et al. 2001); it is also not 
negligible, however, so the question arises as to how the 
two scattering processes should add. Such a calculation, 
while clearly important, is beyond the scope of this pa- 
per. We instead take a pragmatic approach and use two 
different prescriptions: 



Aire p 



max [ £ Z^YyAe,;, (^Ae,imp 



and 



47re^ 



e,ion 



QA, 



e,imp J 



(15) 



(16) 



Here A^ jon is the Coulomb logarithm for a single ion 
of charge number (Z) and mass number {A). Both ap- 
proaches give comparable results in the Debye screening 
limit {9 — > 0) and in the limit {9 3> 1) where impurity scat- 
tering dominates. In the regime where both impurity and 
phonon scattering are comparable, the second prescrip- 
tion (eq. [16]) gives a larger t^^^^ and hence a smaller K. 
For purposes of comparison (§ 4) we compare the conduc- 
tivity of a pure state, e.g., ^°*Ru with that obtained using 
equation (16); this gives the largest variation in K. 

3.4. Sensitivity 

The choice of input physics can dramatically affect the 
relation Teff = 'reff(Tb). The greatest uncertainty lies with 
the calculation of conductive opacities around the melt- 
ing point, F ss Fjvi, and in the crystalline phase. This is 
partly due to our ignorance of the exact composition of 
the envelope. We consider in §§ 4 and 5 different possi- 
bilities for the composition of the heavy elements: ^^Fe, 
^**Ru, and the rp-process ashes. In this section we con- 
sider how a different prescription for the electron-ion con- 
ductivity (Flowers & Itoh 1981; Itoh et al. 1983; Itoh & Ko- 
hyama 1993) than ouj adopted formulae (Potekhin et al. 
1999) would change our results. By looking at the sen- 
sitivities of our results to the choice of conductivity, we 
can understand how our results vary in response to the 
general uncertainties in input physics. 

Figure 2 highlights this problem. Away from the melt- 
ing point in the liquid regime, the fitting formulae given 
by Potekhin et al. (1999), Kpot., and Itoh & Kohyama 
(1993), Kitoh/ are in good agreement, |1 - Kp„t./Kit„h| « 
5 %-60 %. The agreement begins to fall apart near the 
melting point and the conductivities in the crystalline 
regime differ by a factor of 2-3 (for a discussion, see 
Potekhin et al. 1999). Because our conductivity is uncer- 
tain, in any case, for a multi-species plasma (§ 3.3), we 
consider how variations in the thermal conductivity af- 
fect the relation between Ty, and Teg. 
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For a neutron star envelope of fixed composition with a 
given Teff, the temperature profile T(y) can roughly be di- 
vided into three regions: the radiative zone, Kc :§> k^; the 
sensitivity strip, Kc ~ k^; and the isothermal zone Kc » ftr- 
The sensitivity strip is so named because changes to the 
conductive opacity in this region strongly affect the tem- 
perature profile of the envelope (Gudmundsson et al. 
1983; Potekhin et al. 1997; Ventura & Potekhin 2002). A 
change in the conductive opacity in the sensitivity strip 
changes the region where the sensitivity strip lies. Since 
this region controls the transition from a power law ra- 
diative solution to a isothermal zone solution, it is critical 
that the conductive opacity be well-understood here. 

The location of the sensitivity strip yss is roughly where 
Kc ~ Kr. Setting Kc — Ck^ (where C is an arbitrary con- 
stant that contains our uncertainty regarding the compo- 
sition and scattering integrals in Kc and Kr), using an ideal 
gas equation of state, and setting factors of order imity to 
unity, we find that 

Q-l/3 jl7/6 



! 2.2 X 10^ g cm" 



(17) 



(Z + l)gl4 ' 

To relate i/ss to T^ff, we insert the solution to the flux equa- 
tion (eq. [1]) in the radiative zone. Since the dominant 
opacity is from free-free absorption, we can take as our 
opacity k oc pT^'^l'^. Inserting this into equation (1) and 
again using an ideal gas equation of state, we find that 
oc T^^/'^. Solving for T and inserting all the appropriate 
numerical factors, we have 



0.16y' 



,4/17 



A(Z + 1) 



2/17 



(18) 



hiserting the expression for yss, equation (17), into equa- 
tion (18), one finds that the temperature in the sensitiv- 
ity strip scales as Tss oc C"*/^^. To estimate how Teff scales 
with the microphysical input, we note that if Tss « Tt, then 
Teff oc C^/^. Therefore, this calculation is moderately sen- 
sitive to fractional xmcertainties of order 10 % in the input 
physics. Since the sensitivity strip is in the regime where 
the envelope matter remains a liquid (Ventura & Potekhin 
2002), our sensitivity to the microphysics of the crystalline 
matter is relatively small, especially for the surface tem- 
peratures (Teff ~ 10^ K) of interest. 

4. THE VARIATION OF EFFECTIVE TEMPERATURE 

Having laid out our microphysical tools, we are now 

ready to explore how the changing envelope stratifica- 
tion varies the relation between the deep crust temper- 
ature Tb and the effective temperature Teff. To do this, we 
adopt a two-layer model with a variable column depth 
yi of the top layer. The outer layer is composed of H 
or "^He, and the inner layer either pure ^''Fe, ^""^Ru, or 
ashes from rp-process burning. We integrate equation (1) 
numerically using an Adams predictor-corrector method 
(Hindmarsh 1983). As a boundary condition for equa- 
tion (1), we apply the Eddington approximation^ at the 
photosphere, k{p, Teff)yph — 2/3. For a given Teff, we then 
integrate equation (1) inwards to yb — 10^^ g cm"^. At this 
column, the thermal gradient becomes nearly isothermal, 

^ The value of for a given Tgff is insensitive to the precise location of the 



and Tb = T(y — yb) is approximately the interior temper- 
ature. The inverse relation Teff(Tb) is then found by itera- 
tion. As a check, we compared ouj calculations to those 
of Potekhin et al. (1997) for a ^^Fe envelope and a "fully 
accreted" envelope (H/He/C/Fe layers). For a given Tb, 
the fractional difference between our value of Teff and 
that computed from the fitting formula of Potekhin et al. 
(1997) is of order 5 %, with the largest deviation occurring 
when the H/^''Fe interface is in the sensitivity strip. 

To illustrate how the opacity changes with the variation 
in the location of the interface, we show in Figure 3 a two- 
layer neutron star envelope (H superincumbent on ^"^Ru) 
with yi = lO'* '* g cm"^ (panels a, c) or y\ = V^-^^ g cm"^ 
(panels b, d). In both cases Tb = 7.5 x 10^ K. The two 
top panels (a, b) depict the temperature, while the bottom 
panels (c, d) display the total opacity (solid line), radiative 
opacity (dotted line), and conductive opacity (daslied line). 
When the interface is at a low column, both the radiative 
and conductive opacities play a role. At higher column, 
the conductive opacity dominates at the location of the in- 
terface. In both cases there is a substantial increase in the 
opacity of a ^°*Ru layer from a H layer, reflecting the in- 
crease in the ion charge for bremsstrahlung and electron- 
ion scattering. 

At low densities, the opacity is dominated by radiative 
processes (mostly free-free). For a free-free dominated en- 
velope, T(y) oc y'^l^'^ (eq. [18]); as a result, along the trajec- 
tory {y, T(y)}, the free-free opacity is oc y^^^^^ (Fig. 3, 
dotted lines) and is nearly constant. As the composition 
changes from H to ^°*Ru, the opacity jumps by a factor 
w 44^/104. Where the electrons are degenerate and the 
heat transport set by electron conduction, dT/dy — > 0, 

and the Gaunt factor scales as n^^^'^, so kh oc y^/*. In the 
limit where the electrons are degenerate and relativistic, 
the electron conductive opacity scales as Kc oc [T(y)]^/ y, 
and T(y) is nearly constant. The jump at the interface is 
Akc < 7? I A = 4A^/10A because the stronger ion-ion corre- 
lations (parameterized by F) decrease the scattering rate 
T^J and offset the increase from the larger ionic charge. 

Figure 4 shows the emergent flux (T^g) for a two-layer 
neutron star envelope as a fimction of yi. The top layer 
is either pure H (thin lines) or ^He {thick lines), and the 
bottom layer is composed of either ^^Fe (solid lines), ^°*Ru 
(dotted lines), or an rp-process mixture with the mean-ion 
approximation (eq. [16]; dashed lines). For each composi- 
tion, we used 4 bottom temperatures: Tb = 3.75 x 10^ K, 
7.5 x 10^ K, 1.5 X 108 and 3.0 x 10» K. The difference 
in Tgg between an outer H layer and an outer *He layer is 
t5^ically smaller than the difference between the differ- 
ent ash compositions and locations of the interface. Were 
we to calculate Te^on using equation (15) rather than equa- 
tion (16), the resulting curves would have lain between 
those of an ^"^Ru layer (dotted lines) and the rp-ashes with 
the mean-ion approximation (dashed lines). The emer- 
gent flux is insensitive to the location of the interface for 
y £ 10^ g cm"^; as the interface moves deeper the lower 
opacity of H layer reduces |dT/dy| so that Teff increases 
for a fixed Tb. Note also that the profile for the multi- 
photosphere, so this approximation is sufficient for our purposes. 
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species ash has a lower Tgff (higher opacity) than that of 
a pure ^"''Ru layer despite having a smaller (Z). Because 
the Fy for each species in a multi-component mixture is 
reduced relative to that for a pure species, A^j is larger at 
high densities, where <1. 

By how much can the depth of the outermost layer 
vary? The maximum depth of the H layer is set by the 
reaction p(e^,i^)n. This reaction occurs for e-p > nin — 
mp+me = 1.29 MeV, or p w 1.3 x 10^ g cm"^. With our 
assumed surface gravity, the corresponding colurrm is 
y — 1.6 X 10^° g cm^^. For a pure *He layer, the maxi- 
mum depth would be where the strongly screened 3a re- 
action (Fushiki & Lamb 1987) ignites: p > 10* g cm^^ for 
temperatures < 10* K. Accretion to this depth requires 
a very slow accretion rate over a long time. For the col- 
umn accretion rate needed to power the quiescent ther- 
mal emission, m w 1 g cm^^ s^^, the time needed to ac- 
crete H to the electron capture depth is 500 yr. This could 
possibly occur for long-recurrence time transients such 
as KS 1731-260 (Rutledge et al. 2001c; Wijnands et al. 
2001). For short recurrence-time transients, such as Aql X- 
1, the H layer cannot be appreciably thicker than where 
*He unstably ignites, i/ign w 10* g cm^^ (see Bildsten 1998, 
and references therein). How thin the light-element layer 
might be is more difficult to determine. As noted in §2, 
the column of H deposited after the last burst for spheri- 
cally symmetric accretion is < i/ign. If the accretion were 
to occur onto only a small fraction of the surface and then 
later spread over the surface, than the residual coliram 
could be much less than i/ign . 

5. AQL X-1, CEN X-4, AND SAX J1808.4— 3658 

Having explored the variation induced in Teg by vary- 
ing i/i and the composition of the envelope, we now de- 
scribe the thermal structure of Aql X-1, Cen X-4, and 
SAX J1808.4-3658. Aql X-1 and Cen X-4, despite hav- 
ing similar binary orbital periods (19 hr and 15.1 hr, re- 
spectively), have very different outburst morphologies: 
Aql X-1 goes into an w 30 d outburst on a roughly yearly 
basis, while Cen X-4 has had just two recorded outbursts 
(only one of which contributed significantly to the total 
observed fluence) in the past 33 yr. SAX J1808.4-3658 is 
distinguished from both Cen X-4 and Aql X-1 by virtue 
of having pulsations (Wijnands & van der Klis 1998) in 
the persistent emission; its orbital period is also much 
shorter (2.01 hr; Chakrabarty & Morgan 1998), and it 
quite possibly accretes from a substellar mass companion 
(this also explains its low time-averaged accretion rate, 
(M) - 10"" M© yr-i; Bildsten & Chakrabarty 2001). Fig- 
ure 5 displays a summary of calculations for these three 
objects: from top to bottom, Aql X-1 (T^u = 2.1 x 10*' K), 
Cen X-4 {T^s = 8.8 x 10^ K), and SAX J1808.4-3658 (T^s = 
6.8 X 10^ K). We fix the composition of the outer layer to 
be ^He (solid line) with a column yi — 10* g cm^^, and vary 
the composition of the inner layer between ^^Fe [dotted 
line), '^^^Ru {dashed line), and rp-process ashes {dot-dashed 
line). We now explain each calculation in more detail. 

Chandra observations of Aql X-1 (Rutledge et al. 2001a) 
find that k^T^g ^ = 135^]^* eV. For a fiducial neutron star 
of mass M = 1.4 M© and R = 10 km, the proper effec- 
tive temperature is Teff = Teff,oo(l + z) = 2.1 x 10^ K. At 



this temperature and y^, the fractional difference in Tb be- 
tween an outer layer composed of H and one composed 
of ^He is rb(H)/rb(^He) - 1< 0.02. Using Figure 4, the 
variation in the emergent flux induced by a changing yi 
is, at this Tb, roughly a factor of 2; changing the composi- 
tion of the ash from ^''Fe to rp-process ashes increases T^g 
by a factor w 1.6. 

This calculation is consistent with widely spaced ob- 
servations of Aql X-1, which find variability by a factor 
of 1.9 over a timescale of years (i.e., over several out- 
burst/recurrence timescales). We cannot explain changes 
on timescales less than this, however. As mentioned in 
§ 1, Rutledge et al. (2002a) observed Aql X-1 with Chan- 
dra on four successive occasions following the Novem- 
ber 2000 outburst. The intensity was observed to de- 
crease by a factor of w 0.5 over 3 months and then in- 
crease by a factor of w 1.4 over 1 month. In addition, 
short timescale (< 10^ s) variability was evident in the last 
observation. Our present calculation does not address 
the short timescale variability. The fact that the intensity 
first decreased and then increased also cannot be accomo- 
dated in the scenario we outline in this paper. Moreover, 
the timescale for sedimentation (eq. [9]) is far too brief. 

For Cen X-4, Chandra observations find that k^Teff — 
(76 ± 7) eV (Rutledge et al. 2001b). The last known accre- 
tion outburst occurred in 1979. Assuming no outbursts 
have occurred since then, the crust has had several ther- 
mal times to relax. As with Aql X-1, we again set the 
outer layer to be pure ^He (Fig. 5; solid line). The reason 
we choose the outer layer to be *He reflects our preju- 
dice that H is consumed as the accretion rate decreases 
at the end of the outburst. Changing the outer layer 
to H increases Tb by a factor rb(H)/rb(*He) - 1 < 0.08. 
For transients with recurrence times of decades or longer, 
such as Cen X-4 and KS 1731-260, it is possible that 
residual accretion could substantially increase the depth 
of the H/^'He layer. Wijnands et al. (2001) and Rutledge 
et al. (2001c) foxmd the luminosity of KS 1731-260 to be 
w 3 X 10^^ ergs s~^; this constrains the quiescent accretion 
rate to < 2 x 10"^^ M© yr"^. Accretion at this limiting 
rate increases yi by 10* g cm^^ every 30 years. By itself, 
this can change the brightness by a factor of 1.2 on this 
timescale. This effect may be dwarfed by the thermal re- 
laxation of the crust, however, which also occurs over a 
timescale of decades (Rutledge et al. 2001c). A simula- 
tion of the time-dependent thermal luminosity for such 
sources is beyond the scope of this introductory paper; 
for now we just note this interesting possibility. 

Finally, we turn our attention to SAX J1808.4-3658. 
This source is rather dim in quiescence (Stella et al. 2000; 
Dotani et al. 2000; Wijnands et al. 2002). We estimate the 
surface effective temperature from the flux reported by 
Wijnands et al. (2002); note that this is not a bolometric 
flux, and may also include a contribution from a power- 
law. For the estimated temperature, the overall stratifica- 
tion is similar to that of Cen X-4, with Tb w 2.8 x 10^ K. 
The variation in T^a is approximately a factor of 4 at 
Tgff < 10^ K. Given its short recurrence time and cold Teff, 
SAX J1808.4-3658 may be the best source from which to 
observe the effect described in this paper. If the previous 
few years are typical, then either Chandra or XMM can 
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likely observe different quiescent epochs. 

6. SUMMARY AND DISCUSSION 

In summary, we find that (1) the timescale for the neu- 
tron star envelope to segregate into layers is much less 
than the outburst recurrence time, making the surface ef- 
fective temperature sensitive to the mass of H/^He re- 
maining on the surface at the end of an outburst; (2) vari- 
ations in the composition of the heat-blanketing envelope 
can lead to variability, by a factor of 2 to 3, in the ther- 
mal quiescent flux from neutron star SXTs; (3) the crust 
temperatures of Aql X-1, Cen X-4, and SAX J1808.4-3658, 
for the quoted effective temperatures, are respectively 
3.3 X 10^ K, 4.9 X 10^ K, and 2.8 x W K wifli a fractional 
uncertainty, from yi and the ash composition, of roughly 
20 %. 

The measured T^g and inferred of Aql X-1 have in- 
teresting implications for the interior temperature of the 
neutron star. As described in Brown (2000), when the 
temperature in the crust is sufficiently hot, there is an 
inversion of dT/dr: heat flows inward from the crust to 
the core, where it is balanced by the neutrino luminosity 
L^. A calculation similar to that described in Brown & 
Ushomirsky (2000), with the "moderate superfluid" case 
discussed in Brown (2000) finds that the ratio of quiescent 
photon luminosity to neutrino luminosity is Lq/L,, = 3.6. 
Our choice of superfluid transition temperatures in these 
calculations are similar to those employed by Yakovlev 
et al. (2001) in fitting to observations of cooling neutron 
stars. This suggests that for sources with higher (L) than 



Aql X-1, the simple distance-independent relation be- 
tween the quiescent flux and the time-averaged outburst 
flux (Brown et al. 1998) does not exactly hold. For sources 
with lower (L), neutrino cooling is not important, imless 
an enhanced (beyond modified Urea; see Colpi et al. 2001) 
neutrino emissivity operates. A direction for future work 
would be to incorporate recent fits (Kaminker et al. 2001) 
of cooling neutron star observations. 

We reiterate that the best place to observe the effect of 
a changing light-element layer mass is a source such as 
SAX J1808.4-3658, which has both a low T^a (Wijnands 
et al. 2002) and a short recurrence timescale (2 outbursts 
in the past 6 years). At lower Tb, the effective temperature 
can vary by as much as a factor of 4, while the frequent 
outbursts allow comparison between different envelope 
layer ings for a fixed core temperature. 
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Fig. 1. — Schematic of different physical regimes for a quiescent neutron star envelope. We show the conditions F = 1 and F = Fm (solid lines), 9 = 1 
{dotted lines) and t/j = 10 (dashed lines), for both H (top panel) and '"^Ru (bottom panel). Circles denote the envelope structure for a two-layer (H over 
^^^Ru) envelope at T^g = 2.1 x 10^ K (top curve) and T^n = 8.8 x 10^ K (bottom curve), as appropriate for Aql X-1 and Cen X-4, respectively (see text). 
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Fig. 2. — Ratio of the conductive opacity as fitted by Potekhin et al. (1999), Kpot./ to thiat fitted by Itoh et al. (1983), Knoh/ ft"' various plasma 
compositions as a function of density at a temperature T = 10* K. The discontinuities occur at the liquid-solid phase transition. They reflect our 
choice to use (for the results of Itoh et al. 1983) the liquid metal approximation for r < r^j and the electron-phonon scattering for r > Ty^. 
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Fig. 3. — Temperature and opacities for a two-layer (H over ^"^Ru) envelope. The interface between the layers is denoted by the thin solid vertical 
lines: i/j = 2.5 x 10'' g cm^^ {left panels) and i/j = 1.3 x 10* g cm^^ {right panels). The top panels depict the temperature; the bottom panels depict the 
total opacity {solid lines), radiative opacity {dotted lines), and conductive opacity {dashed lines). 
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Fig. 4. — Variation of T*jj with interfacial column i/; for a two-layer envelope with a fixed inner boundary temperature Tj,. The top layer is either H 
(thin lines) or *He (thick lines); the bottom layer is either ^^Fe (solid lines), ^"^Ra (dotted lines), or a very impure (Q = 233; K calculated in the mean-ion 
approximation) mixture from an rp-process burst (dashed lines). Each group of curves corresponds to T^, = 3.75 x 10^ K, 7.5 x 10^ K, 1.5 x 10« K, and 
3.0 X 10» K. 
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Fig. 5. — Thermal structure of a two-layer neutron star envelope with the interface between the layers at yj = 10^ g cm ^. The three structures, from 
top to bottom, are for Teff = 2.1 x 10^ K (Aql X-1), Teg = 8.8 x 10^ K (Cen X-4), and = 6.8 x 10^ K (SAX J1808.4-3658). The composition of the 
outer layer is *He {solid line). The composition of the inner layer is either ^^Fe (dotted line), ^"^Ru {dashed line), or rp-process ashes with the mean-ion 
approximation (eq. [16]) for the conductivity {dot-dashed line). 



